Absence of the diffusion pole in the Anderson insulator 
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We discuss conditions for the existence of the diffusion pole and its consequences in disordered 
noninteracting electron systems. Using only nonperturbative and exact arguments we find against 
expectations that the diffusion pole can exist only in the diffusive (metallic) regime. We demonstrate 
that the diffusion pole in the Anderson localization phase would lead to nonexistence of the self- 
energy and hence to a physically inconsistent picture. The way how to consistently treat and 
understand the Anderson localization transition with vanishing of the diffusion pole is presented. 
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The Anderson insulator is a specific disordered or 
amorphous solid with available quasiparticle states at 
the Fermi surface but with no diffusion or charge trans- 
port at long distances. Possibility of the absence of dif- 
fusion in impure metals and alloys was proposed by P. 
W. Anderson on a simple tight-binding model of disor- 
dered noninteracting electrons Jj. Since then, vanishing 
of diffusion, now called Anderson localization, has at- 
tracted much attention of both theorists and experimen- 
talists 01 • In spite of a considerable portion of amassed 
experimental data, disclosed various specific and general 
aspects of the Anderson metal-insulator transition and 
a number of theoretical and computational approaches so 
far developed, we have not yet reached complete under- 
standing of Anderson localization. Although the basic as- 
pects of the critical behavior at the Anderson localization 
transition are known, the position of this disorder-driven 
metal-insulator transition within the standard classifica- 
tion scheme of phase transitions with control and order 
parameters remains yet unclear. 

In this Letter we reexamine the description of Ander- 
son localization with averaged functions. Configurational 
averaging is an important means to restore translational 
invariance in disordered systems and to provide repro- 
ducibility (sample independence) of the findings. Disor- 
dered systems after averaging behave as pure ones with 
effective correlations between the motion of individual 
quasiparticles. The translationally invariant description 
is hence the proper tool for comprehending the Anderson 
localization transition in the way we understand phase 
transitions in clean interacting systems. 

The aim of the Letter is to set forth general con- 
straints on the translationally invariant description of 
the Anderson localization transition resulting from ex- 
act relations, conservation laws and equations of motion 
for the averaged Green functions. We explicitly demon- 
strate that the weight of the diffusion pole cannot be fixed 
to unity as dictated by the Ward identity and that the 
diffusion pole must vanish in the localized phase. The 
diffusion pole in the localized phase would lead to nonex- 
istence of the self-energy and to an unsolvable equation 



of motion for the two-particle vertex. 

The fundamental building blocks of the transla- 
tional invariant description of disordered electrons are 
averaged one and two-particle resolvents G(k, z) and 
^kk' (-^+ ' ; ' respectively. Here k, q are fermionic 
and bosonic momenta and z± are complex energies. We 
use the electron-hole representation for the two-particle 
Green function with k and k' for incoming and outgoing 
electron momenta. The bosonic momentum q measures 
the difference between the incoming momenta of the elec- 
tron and the hole. 

Fundamental functions for the description of a linear 
response of an electron gas to an external electromagnetic 
perturbation are the density response and the electron- 
hole correlation function defined from the two-particle 
resolvent as 

X(q, iVm) = X! X! ^kk' + q), 

kk' n= — ca 
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For convenience, the density response x defined at 
the temperature axis, while the electron-hole correlation 
function for real frequencies. The superscript in- 

dices i?, A relate to the limits of complex energies from 
which the first and second real energy variables were 
reached. 

It is useful to go over from the full two-particle resol- 
vent to a two-particle vertex F defined as 



'-'kk'l^^ 



,z_;q) = G(k-f-q, z_ 
-rkk'(^+,2~;q)G(k' 



)G(k,z_)[l 

Kq,z+)G(k',z_)] . (2) 



The two-particle vertex obeys a Bethe-Salpeter equation 
of motion. Since we use the electron-hole representation, 
the most natural way to construct the Bethe-Salpeter 
equation is to use ladders in the electron-hole scattering 
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channel. We then can write 



rkk'(q) 

1 



N ^ 

k" 



kk'(q) 

A^^(q)G+(k" 



q)G_(k")rk''k'(q) (3) 



where we introduced the electron-hole irreducible ver- 
tex A'^''. We suppressed the energy variables in Eq. (Q, 
since they are not dynamical quantities and act only as 
external parameters. They are easily deducible from the 
one-electron propagators G±(k) = G(k, z±). 

We have up to now used only definitions or represen- 
tations of the two-particle functions. The first impor- 
tant equation demanding a proof is the VoUhardt-Wolfle 
identity expressing conservation of the norm of the wave 
function of a free electron scattered on a static impurity 
potential 

k' 

X [G^(k',S + co)-G^(k',S)] . (4) 

Here we denoted T,^, T,^ the retarded and advanced self- 
energy and K^^,{E+uj, E) = M^^,{E+uj+iQ+ , E~iQ+;^. 
Equation (0J was proved diagrammatically in Ref. 
A nonperturbative proof exists trough equivalence of 
Eq. 10} with the Velicky identity 0, . We hence can 
consider Eq. 10} as a nonperturbative property of solu- 
tions of the Anderson model of disordered electrons . 

The first fundamental consequence of Eq. Q is the 
existence of the diffusion pole in the electron-hole cor- 
relation function <1>^^. Introducing a quantum diffusion 
function D(q, lj) we can represent the density response 
in the hydrodynamic limit (jf ^ as 



x(q,w) 



—iuj + D{Lj)q'^ 



(5) 



where xo = limq^o limt^^o x(q, t^^) and D{uj) = £'(0,a;). 
The dynamical diffusion coefficient D{ljj) reduces in the 
static limit to the diffusion constant related at zero tem- 
perature to the static optical conductivity by the Einstein 
formula D = D(0) = a/e^np-, where n^- is the density of 
states at the Fermi energy 

Using Eqs. |^ we can now expand the density re- 
sponse in the hydrodynamic {q 0) and quasistatic 



(w — > 0) limits maintaining D{ui)q 



We find that 



its leading term at zero temperature is governed by the 
electron-hole correlation function 

X(q, u) = x(q, 0) -f ^ ($i^(q, u) + Oiq'^)) + 0{u) . 

(6) 

Inserting the asymptotic expansion Eq. ^ in Eq. |(S} we 
derive an explicit representation of the diffusion pole as 



the leading low-energy term in the electron-hole correla- 
tion function 



*i-'(q,^) = 



2'Knp 



-iLU + D{uj)q'^ 



+ 0{q\^°). (7) 



Although the denominator in the density response in 
Eq. (0 is identical with that from the electron-hole cor- 
relation function in Eq. |7}, it is only the latter function 
that is singular in the low-energy limit. This singularity, 
the diffusion pole, is a consequence of the particle-number 
conservation expressed mathematically in the Ward iden- 
tity 0}. The only parameter depending on the charac- 
ter of the underlying solution in Eq. (|7} is the dynami- 
cal diffusion coefficient D{uj). This simplified low-energy 
asymptotics of the electron-hole correlation function is 
the basis for the existing scaling and field-theoretic ap- 
proaches to the problem of Anderson localization Q . 

According to the above reasoning the diffusion pole 
must exist also in the state with localized electrons, where 
the dynamical diffusion coefficient is expected to vanish 
in the static limit as yl 



D{uj) 



(8) 



where ^ is the localization length. We, however, show 
that such a picture is not sustainable, since the diffusion 
pole in the localized phase would hinder the existence of 
a finite self-energy almost everywhere within the energy 
band. 

To manifest this we first show that in systems invari- 
ant with respect to time reversal the diffusion pole enters 
the electron-hole irreducible vertex A'^^ as the so-called 
Cooper pole. It displays the same low-energy asymp- 
totics as the electron-hole correlation function <I>^^. 

Time inversion in nonintcracting systems with elas- 
tic scatterings only is expressed as inversion of the 
particle propagation k — > — k, i. e., the electron 
and hole interchange their roles. The time-reversal 
invariance then means that G(k, z) — G(— k, z). 
For the two-particle resolvent we then obtain either 
G'kk'(^+' q) = G%-'kiz+,z-;q + k + k') or equiva- 

(2) (2) 

lently Gj^^,(z+,z_;q) = Gj^^,(z+, z_; -q-k-k') leading 
to Fkk'Cq) ~ r_k'-k(q + k -I- k'). However, the electron- 
hole irreducible vertex A^^ is not invariant in the same 
way the full two-particle vertex F is. We actually have 
Ai'j^,_i^(q + k + k') = Ag.|.,(q), where A''" is the electron- 
electron irreducible vertex being generally different from 
the electron-hole one. The Bethe-Salpeter equation ||3} 
after time reversal is hence transformed to another Bethe- 
Salpeter equation with the electron-electron irreducible 
vertex and modified momentum convolutions. The ex- 
istence of several Bethe-Salpeter equations for the same 
full two-particle vertex F is a well known ambiguity in 
the definition of the two-particle irreducibility |5j . 

The low-energy singularity in the electron-hole cor- 
relation function is a consequence of the existence of a 
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singularity in the full two-particle vertex T^^, since the 
former function is the latter one multiplied with one- 
electron propagators and summed over fermionic mo- 
menta, cf. Eqs. ((rE|) and It means that r^j^(q) 
diverges with w, — > for almost all fermionic mo- 
menta k, k' so that its singularity survives in as 
the diffusion pole, Eq. However, due to the electron- 
hole symmetry, the full two-particle vertex must show the 
same singularity also in the limit ui, (k + k' + q)^ 0. 

Further on, the (leading) singularities in the full ver- 
tex r are already contained in the irreducible vertices 
A*^'' and A'^'^. This conclusion follows from the so-called 
parquet equation. It is an expression of topological 
nonequivalence of different two-particle irreducibilities. 
That is, the reducible function of one type (convolution 
of two or more irreducible vertices of the same type) is 
irreducible in the other irreducibility channels. In case of 
two irreducibility channels (two types of Bethe-Salpeter 
equations) we can write the parquet equation as 



rkk'(q) = Ark'(q) + A£'^,(q) - Xkk'(q) 



(9) 



where J = A^'* n A^'^ is a vertex irreducible in both irre- 
ducibility channels (completely irreducible vertex) ||^. 

We deduce the form of the pole in the electron-hole 
irreducible vertex A'^'* from the asymptotic limit to high 
spatial dimensions 9]. There we obtain that the com- 
pletely irreducible vertex 2 is local and regular. The 
singularity of the full vertex in the limit w,q^ -^0 
is contained in the electron-electron irreducible vertex, 
A^*^, while the singularity for w, (k -f- k' -I- q)^ ap- 
pears only in A'^^. The low-energy singularities of the 
irreducible vertices are fixed in momentum space by the 
asymptotics d — > oo if this limit is analytic for the two- 
particle functions. Since the neglected contributions to 
the two-particle vertex are less singular than the lead- 
ing divergent d — > oo terms containing the diffusion and 
Cooper poles, we can assume analyticity of the high- 
dimensional limit, even if we cannot prove it rigorously. 
Using the high dimensional separation of singularities of 
the full vertex into the irreducible ones we obtain from 
Eq. at zero temperature for q ^ 



RAi 



2TmpX 



-iuj + D{Lu)(k + k')'^ 
(10) 

where A is a measure of the disorder strength. The func- 
tion A^j^ is nonsingular or at least less singular than the 
Cooper pole singled out from the electron-hole irreducible 
vertex, second term on the r.h.s. of Eq. H1U|) . Notice that 
the asymptotic form of the Copper pole in Eq. ifTU)! is ex- 
act only if the the Ward identity is obeyed. The 
Cooper pole in Eq. Hl()|l is the generally anticipated sin- 
gularity of the electron-hole irreducible vertex 

To test consistence of such singular behavior of the 
electron-hole irreducible vertex we define a difference of 



two retarded self-energies 
^We{co) = ^^[E^(k,^-c^)-E^(k,^ + cc;)] . (Ua) 



We will be interested in its low- frequency behavior, i. e., 
the asymptotic limit lim^^o AWsito)- If the self-energy 
I]^{E) is an analytic function of the energy variable E, 
which normally is the case within the energy bands, then 
AWe{oj) must be analytic at w = for almost all en- 
ergies E for which 3E^(i?) > 0. If the self-energy 
is a continuous function of the energy argument, then 
AWe{0) = 0. Finiteness of the derivative of the self- 
energy demands limi^_o |AM/£;((jj)/ti>| < oo, etc. 

We apply the Ward identity Q to represent the r.h.s. 
of Eq. (|lla|l via the electron-hole irreducible vertex. We 
obtain 



N2 



kk' 



[AS^,iE + u;,E) 



ill. 



A^^,{E- 



- LO. 



^^,{E~u,E)\ 3G^(k',i?) 
-uo^E) [G^,{E + uj)-G^,{E)] 
E) [GUE~u:)-GUE)\] 



(lib) 



We use representation (|10|l for the electron-hole irre- 
ducible vertex to derive the low-frequency asymptotics 
lim^^^^o AW£;(w). Actually we are interested only in the 
most singular part being generated by the Cooper pole. 
It is a straightforward task to perform momentum inte- 
gration on a d-dimensional hypercubic lattice leading to 

m 









D{u)kl 


u 






D{u)kl 



d/2 
d/2 



In 



D{u. 



for d ^ 4Z, 



for d = Al. 



(12) 



In the diffusive phase, D{uj) = D > 0, that reahzes 
only for d > 2, we find nonexistence of derivatives of the 
self-energy independently of the energy position within 
the band. Realization of a self-energy with nonexisting 
(diverging) [{d — 2)/2]th derivative everywhere is appar- 
ently unphysical for d < oo. 

The situation worsens if we evaluate limi^^o AWEit-^) 
for the localized solution with the low-frequency dynami- 
cal diffusion function from Eq. © . We then obtain inde- 
pendently of the spatial dimension AW^^^ {w) ^ Cd/oj. 
This is a catastrophe indicating nonexistence of the re- 
tarded self-energy as a well defined mathematical func- 
tion. This may not happen. Hence, our picture based on 
the Ward identity Q resulting in the Cooper pole in the 
electron- hole vertex, Eq. (|HJ|I . is inconsistent. The real 
solution of the Anderson model must behave differently 
from the above derived singular behavior. 
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The correct low-energy behavior of the electron-hole 
vertex can be assessed from the exact asymptotic solu- 
tion in high spatial dimensions. The limit to high spatial 
dimensions leads to simplifications in momentum convo- 
lutions so that one can solve the parquet equations for the 
irreducible vertices K'^^ and A^*^ asymptotically exactly 
0. In this asymptotic solution the Ward identity 1^ 
is not obeyed for nonzero real energy differences uj ^ Q. 
A deviation from the Ward identity leads to a modified 
diffusion (Cooper) pole. In the localized phase we then 
obtain for lo, (k -I- k')^ — *■ the electron-hole irreducible 
vertex in the asymptotic form 

where = iD{uj) / A{uo)u) is the square of the localization 
length. The limit limi^^o A{lij)lj becomes finite inside the 
localized phase as well as the localization length ^. We 
can see from Eq. that the electron-hole irreducible 
vertex does not display a low-energy singularity. It means 
that there is also no diffusion pole in the electron-hole 
correlation function in the localized phase. This is a nec- 
essary condition for the existence of the self-energy as 
a well behaved function along the real axis of the energy 
variable (Fermi energy). 

The fraction |l/A(u;)| < 1 generally determines the 
weight of the diffusion pole. The static value A(0) is fi- 
nite in the diffusive phase, increases with the disorder 
strength and diverges at the Anderson localization tran- 
sition. The renormalized density of states np/AiQ) rep- 
resents the number of diffusive states, that is, extended 
states being able to carry electron charge to long dis- 
tances. These are the only states contributing to charge 
transport, whereas n^^- is the density of all available states 
at the Fermi energy Ep. 

Nonconservation of the averaged particle density in- 
curred by a violation of the Ward identity l@J in the high- 
dimensional solution for the vertex function is a conse- 
quence of the configurational averaging ■ Validity of 
Eq. Q is conditioned by completeness of Bloch waves for 
all relevant configurations of the random potential. Dur- 
ing the configurational averaging we lose the spatially re- 
stricted localized states in the thermodynamic limit. To 
form a representation Hilbert space we have to fix bound- 
ary conditions for the eigenstates of the random Hamilto- 
nian. The extended states vanish at the volume boundary 
as while the localized states vanish exponentially 

exp{— y/uo}- The two types of states are orthogonal in 
the thermodynamic limit. Restoring translational invari- 
ance restricts our description only to states extended over 
the whole volume. The only localized states surviving 
configurational averaging are the periodically repeated 
ones. If there is a macroscopic portion of configurations 
with spatially bounded states breaking translational in- 



variance near the Fermi energy, the number of available 
translational invariant states is less than the total num- 
ber of states. Hence, the probability to find the electron 
in an extended state is less than one. 

To conclude, we showed that the diffusion pole in 
the electron-hole correlation function cannot survive be- 
yond the mobility edge where all states become localized. 
We used only exact and nonperturbative arguments and 
found incompatibility of the Ward identity imply- 
ing a singular low-energy asymptotics of the electron- 
hole correlation function, Eq. (0, with the existence of 
the retarded self-energy as a well defined function of the 
energy variable. We demonstrated that a necessary con- 
dition for the existence of a regular (finite) self-energy 
in the localized phase is regularity of the electron-hole 
irreducible vertex in the low-energy limit. Supported by 
the asymptotic solution of the Anderson model in high 
spatial dimensions we concluded that the weight of the 
diffusion pole in the metallic regime is 1/A(0) < 1. This 
weight decreases with increasing disorder strength and 
vanishes at the localization transition. The most strik- 
ing consequence of our finding is that it is not the diffu- 
sion constant D{{)) but the weight of the diffusion pole 
1/^(0) < 1 that is the relevant control parameter for the 
Anderson localization transition. 
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